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ABSTRACT 


SPACE AND ENERGY DEPENDENT NEUTRON 
SLOWING DOWN TIME IN GRAPHITE 


Subnitted to the Department of Nuclear Engineering on 
peeuemoer <6, 1967 in partial fulfillment o7mece CC 


for Gne deg 


mee Of master of Science 


Peo LOOu Cube, graphit¢ moderating assembly was con- 


merucevea at the M.1L.t. 


Using the pulsed neutron experimental installation, a neu- 
tron burst of 200 nsec width was injected into the center 
of the assembly. Prompt gamma radiation from a neutron 
Beesonance capture by an indicator foil in the moderator 
mas measured, with a 3 inch x 3 inch sodium iodidé scin-— 


feliation 
tron pulse. 


detector, as a function of time after the neu- 


The energy dependence of slowing down time 


was experimentally measured by using resonance absorber 
foils of cadmium, indium, and gold with isolated reson- 
ances of 0.2, 1.46 and 5.0 eV respectively. By vary- 

ing the distance of the foil from the source, the spatial 
dependence of slowing down time was also investigated. 


Comparisons are made with slowing down theory. 


mental results, together with their estimated uncertainties, 


are tabulated below: 


ec 
12.5 
18.0 
28.0 
ie? 
44.0 
4A 5 


; 


t. (usec) 

0.2 1.46 5.0 
47.50 + 3.00 19. 75 oe eee 
a> eee 19.37 * 225 10.10 + 
47.25 + 2229 19.00 + .50 10.87 + 

- 20.00 + .50 oe 
= a 7S 12.50 + 
HON?) A aeas 21.62 + 50 a 
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CHAPTER 


Introdiction 


In his M.I.T. Ph.D. thesis, Adamantiades (1) demon- 
strated eee pulsed neutron technique using a resonance 
absorber as a point indicator could be used to measure sys- 
ten transient effects such’ as the slowing down time. He 
designed-and installed an experimental installation for 
pulsed meutron experiments, and made some measurements in 
water. Vaughan (2) refined Adamantaides’ system, and measured 
iie Spatial dependence of slowing oun time in water for 
three different resonance absorbers. Moller has carried 
Ou extensive experiments in the field of Jou slowing 
down for both light (3,4) and heavy (5) water. However, 
meme does not appear to be any published measurements of 
meg@enslowing down time for heavier moderators. Therefore, 
Meeoois present experiment, the pulsed technique has been 
fended to investigate slowing down time in graphite. 

In the pulsed neutron method, a short monoenergetic™ 
bes, Of neutrons is injected into the moderating medium 
meme DOint. The neutrons, originally high energy, gradually 
Flow down by elastic collisions with the moderating medium. 
tf a neutron indicator is incorporated into the moderating 
medium, it is possible to measure the time difference between 


pulse injection and neutron arrival at the indicator. A 





time distribution is formed from many time differences, 
measured under wae CXACt Same congue ron. 

Since most materials emit prompt gamma rays upon the 
capture of a Te User Ol, detection of these gamma rays pro- 
mvides one metnod by which the time behavior of neutron 
captures may be Pero after the pulse. If an absorber 
with a strong, isolated neutron capture resonance is incor- 
erated into the moderating system as a point indicator, 
then the reaction rate of the neutron flux with the absorp- 
fe 10n resonance is proportional to the time dependent flux 
at that point, for neutrons with energies equal to the 
resonance energy. That is, the neutron dénsity (Cowart lus 
n(x ,E.,t) “ 

N —- neutron density 

© - distance from source 

E.. - energy of absorber resonance 

t -— time 
can be measured. Thus, variation of the position of the 
eeerber and the use of different resonance absorbers should 
yield the time, space, and energy variation of the neutron 
density following the pulsing of the assembly. 

} By using the M.I.?T. Rockefeller Van de Graaff accelera- 
tor as a source of high energy protons, high energy neutrons 


| were obtained by the Li’(p,n)Be’ 


nuclear reaction. The lithium 
Peareet was located at the center of a large assembly of 
graphite, thus approximating a point source in an infi- 


nite medium. The accelerator beam was pulsed by an external. ... 


% 








pulsing system after the protons had been accelerated. The 
pulse width used in these experiments was 200 nsec. Enough 
time was allowed’ to elapse between pulses so that all tran- 
sient neutron effects from one pulse would die out before 
mne next pulse was injected. Absorber foils (2 sed amir 
maaium, and cold, with isolated resonances located at O.2 eve 
1.46 eV, and 5.0 eV respectively, were used to sample the 
Peutron flux at hese energies. A fe vee reel analyzing 
system was Ee ceronodme ee was capable of indicating the 
time difference between pulse injection and the resonance 
capture of a neutron. The most probable time difference 
trom injection to resonance capture was assumed to be the 
neutron slowing down time for that energy at that positions 
The material in this thesis is organized as’ follows: 
Chapter @ is concerned with slowing down theory, the a 
experimental facilities are described in Chapter 3; the 
performance of the experiments is described in Chapter’4: 


and the results are presented and discussed. 
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CHAPTER 2 


Theory 

._ ; 
2.1 General ¢ ; 
§ 


\ 


_ The analytical treatment of pulsed neutron assem- 


Moliles requires a mathematical solution of the general 


transport equation. If an’ expression representing the neu- 


tron density as a function ot bine, energy, and space, and 
satisfying the Boltzmann transport equation with its source 
ana boundary conditions can be found, the behavior of neu- 
trons in tTheysystem can be completely decce Lucien 

Bintortunately, no exact solution of the Boltzmann 
equation, in all its generality, has ever been found: 
but approximate solutions are available. Of these solutions, 
only-a few may be verified because not all quantities can 
pemmicasured experimentally. 

The quantity most readily measured experimentally 


in the slowing down region is the neutron slowing down 


| time (t,). This is the time at which the neutron density 


-at some energy (E) is maximum and can be found by setting 


tne derivative of N(r,E,t) 
N —- neutron density 
Pe —maistance [rom source 
i = neutron energy 


t — time 








Sos 


with respect to time equal to zero. 
It is also Pees to measure the mean slowing down 
PLM Cmte) which is also defined at the first time moment. 
The time eee are integral quantities defined: 

ae oe > = (aces N(v,t) dt/fN(v,t) dt (2a) 

N(v,t) - neutron density as a function of space 
and energy (velocity) | 

t - time 

Both the time to maximum or the slowing down time, 
and mean Berita. down time may be measured axneene ae as 
afunction of space and neutron energy. 

In this experiment, the source may be considered 
isotropic and a delta function in time, energy, and space. 
mieonoush the neutron source energy does vary an the angie 
joethe laboratory system, this source term ean be neglected! 
because it appears in the denominator of the expressions é 
for the time distributions, and it is much greater ne the 
>_> = auevhich the time distribublion 1S measured on 
example, one Be eoe win for the neutron slowing down flux 
in hydrogen 2 

o(v,t) = S-(vt/y,)° -exp(-vt/y.)+(l-v/v,+2y /v,t)+ 

| 28-v./v-8(v-v,)-exp(-v 9 t/Y,) (221. 2)am 
© —- neutron flux 

.vV -— neutron velocity 


S - source strength 





? | i o> | | = 


Y, — Scattering mean free path 
Sie Deweron SOURCE velocity 
t - time 
fowever, because Vo22v, equation (2.1.2) simpli tocmeerer 
o(v,t) = S-(vt/y,)?- exp(-vt/y .) (2.1 
’ Except for positions closer than three mean free 
_ paths away from the physical boundary of the assembly, the 
spatial boundary conditions are those of an infinite medium. 
ior meee tions closer than three mean paths, leakage must 
be taken into account. 
in the energy region lying above a few tenths of an 
eV, the nuclei of the moderating assembly can be considered 
free and at rest. For lower energies, however, the molecular 
structure of the moderator becomes important. 
“Any absorption which occurs in the moderating assembly 
will effect the time distribution. Those neutrons which 
meee tonuszer to slow down, i.e. have more collisions, to 
energy (v) will have a greater probability of being absorbed. 
Thus, absorption contributes to a decrease in ta eta 
graphite, the absorption cross section is very small, and 


absorption effects may be neglected. 


e.¢ Space Independent Slowing Down Theory 


If, in addition to an infinite medium, we assume a 
uniformly distributed source instead of a spatial delta 


function source, the spatial dependence is removed from the 





ee 


Boltzmann EBouations es ine Caquaplom reduccocs ae. 
“1/v+8e(v,t)/8E = -Le.(v)+e,(v) lo(v,t)+ 
+ Je(v?,t)+e(v’—sv) dv'+S-8(v-v,)-8(t) (2.2.1) 
o(v,t) -— neutron flux as a function of space and 
energy (velocity) 
i Ea macroscopic scattering cross section 
2 macroscopic absorption cross section 
—€ - removal cross section 
S$ ~ source strength (n/cm’ -sec) 
Vo 7 source energy 
6( ) - Dirac delta function 
Unfortunately, for a heavier than water moderator, 


fas very difficult to obtain an expression fom the flux, 


Paaenence, for the slowing down time. Using a method 


me introduced by Marshak (1) the time moments can be obtained 


meoemly. A few of the moments, in terms of the moderator 
mass number, obtained by Adamantiades (2) are given in | 
Table i, Of his thesis. Once tne time moments have been 
determined, Pad functions eee found to fit these 
moments. | 
Marshak and vonDardel (3) suggest using the function: 
p(x) = ee Ee) eens yO pce. 2s) 


p=) ormalization constant 


th 


vt/¥ where y, - scattering mean free path 
(M-m)/(M+m) 


M- mass of moderator 


x 


iG) 





De MASS eof neu bom 

b - a parameter 
which ane right asymptotic dependence for the higher 
positive moments. ones varies the parameter (b) with 
(r) while Marshak treats it as a constant. 

For heavier elements, the time dependent neutron den- 
Sioy Gistribution approaches a narrow Gaussian in shape. 
Peis, for heavier elements, i, should approximately equal 
Bo. 

Adamantiades gives a complete discussion of the time 
monent en in Chapter II of his thesis. He also men- 
tions other methods, and gives a number of references deal- 
ing with the problem of a heavier than water moderator. 
None of these references, however, specifically ie Pesulas 


for graphite assemblies. 


me. ss opace Devendent Slowing Down Theory 


Wnen the neutron source is localized, as in this 
eri ment, the effects of diffusion during moderation must 
be taken into account. The problem becomes much more com- 
“plicated since the transport equations 1S now 2 Lumcliomses 
three, rather than two variables. 

As in the space independent case, Adamantiades in Chap- 
ter II of his thesis discusses the theoretical work which 
has been done in this area. This concludes : solution by 
Dyad’kin and Batalina (4), for the space-energy-time neu- 


tron density distribution assuming no absorption, constant 


—<> —i a 
2 m a EE as ——— — 
> 


ot meee, — eee —_—_ = _ 


+ me mee : 





=) | Bee 


scattering cross section, and isotropic scattering. The 
expression involves three factors which are complicated, 
and require coefficients which are dependent on the medium 
being examined. The authors also give a simplified expres- 
foion for te which is: 
tay = CY,/vm> (2-Bp/Vr° +B" py )- 
(1-20/3+(1-Bu/Ve"+P py] (2.3.2) 


i 


as scattering mean free path 
v —- neutron velocity — i 
mT —- average logarithmic energy decrement 
ww - neutron lethargy 
a —- parameter which depends on the moderating medium 
B —- parameter which depends on the TOG a ete emecearn 


i! 
™ — Gistance from source 


Results are given for sand and water mixtures. 
Claesson (5) has also determined the Sete aopendene 
memyron density distribution for both a point and plane | 
source assuming zero absorption, isotropic scattering, 
and constant mean free path. Although his method can be 
extended to moderators of arbitrary mass, curves are only 
given for the case of water at the indium resonance. These 
curves show however, that past 10 cm the slowing down time 
is strongly fependent on neutron energy at birth. 
Adamantiades also develops a method of calculating 
Gime moments using an iterative solution of the stationary 


Geansport equation. The nth time mgment is expressed 


J & 4 





as 


as 
oe pi Bd jp? 2 ae 
up6) = -S(r,vjn) (2.3.3) 
Ho Man /ye ghar}? = O 2a) 
H= (-pV-e,sJav’Sau? 6 (v?—v, py’) ] (2.3.5) 


P,Vsle= vectors for position, velocity since cue 
Bangi. respectively 


E, - TObal macroscopic absorption cross S-ctioen- 


Any calculational method suitable to solve the steady state 


mam be used to solve for the flux moments. 
- David Diamond, a graduate student at M.1I.T., has 
kindly supplied theoretical time distributions at the gold 


Mia indium resonances for a number of spatial points, 
| | ty 


meee COC eGdiffusion approximation. A distribution for each 
=D) oes 


[meat a distance from the source of 18 cm iS Shown in 


Meezure 2.1 (a) and (b). 
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CH Ab ame 


Experimental Equipment 


4.1 General eee 

The equipment necessary for measuring neutron transient 
effects by the PCa eee method can be grouped into Silex 
Catecories: 

1) The neutron source. 

2) A method for pulsing the source. 

3) The moderating assembly. 

44) A ees oei 1Or Ceveccina. analyzing and storing data. 

SA system to monitor the experimental equipment. 

6) A normalizing system to insure coneerane data. 

‘The general layout of the experimental systen showing 
the major components is given in Pigumer so | , 

The neutron source consisted of a target of eel 
lithium, located approximately in the center of the modera- 
Biane assembly, bombarded by 2.4 MeV pretonus yu. produce the 
ml > n)Be? reaction. The high enérgy protons were Brodie 
Mee van de Graaff electrostatic accelerator. The pulsing 
system electrostatically deflected the proton-beam:off the 
1ithium target, allowing it to strike the target for a short 
Pine . 

Thus, the pulse of fast neutrons produced at the tarses 


was ejected into the center of the moderating assembly. 
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These neutrons were then slowed down by elastic collisions 
aT, the medium, and their time distribution yasemeacumed 
by the detecting and analyzing system. This system con- 
sisted of the conventional detector-preamplitier—-amjli {ier 
array, plus electronic equipment for pulse heigaiee moon 
mination and time analysis combined with a multi-channel 
analyzer to store the data. There was a common link between 
Ghe pulsing and the analyzing system which providedus Zero s 
mime reference. : | 

The monitoring system oe used to control the perfor- 

“ 

mance of the pulsing system, and thus to optimize the quality 
of the neutron pulse. The normalizing system yielded the 
time integrated neutron flux for each experiment perforned. 
This information was necessary to properly carry out back- 
ground subtraction. 


A more detailed description of the experimental equip- 


ment will be given in the sections to follow. 


4.2 The Mechanical Layout | 


A’ plan view of the mechanical layout of the system is 


ee ee ae ee oe 
—~ —— 


Even in Figure 3.2. Only a brief description of the sys- 
tem will be Pirro phere wsinee a CAerovuen., Cevalled descrip- 
tion of the system will be given by Adamantiades (1), who 
actually constructed the systen. 

Since the experimental apparatus. was positioned on 


tne floor below the accelerator, the protons were accelerated 
| 
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downward andi bent into the horizontal bean Cube by eoee. ne 
ing and analyzing magnet. The beam tube was kept undersea 
high vacuum by the diffusion pump of the acceleraicm teed 
by a small Vac Ion pump located close to the target. 

The deflecting plates einen a high voltage pulse 
which deflected the beam. The water-cooled chopping slits 
were mounted after the deflection plates for the purpose 
of intercepting the beam when it was swept off the target. 
~The distance between the slits, the shape of the high vol- 
tage pulse, and the deflecting plates configuration deter- 

Mined the shape of the pulse reaching the whence 
' A remotely controlied pneumatic Sets. valve was 
interposed between the chopping slits and the target. 


fy 
The beam could thus be cut off from the target, enabling 


ee 


_ = cee ee me oe 


work to be done around the assembly. The target tube itself 


was constructed of aluminum and was 24 inches y toncwmes, e 
mom O.D., with a 1/16 inch wall thickness. The target 
assembly was Located at the end of the target tube. 


It was decided to cool the target in order to protect 


moet rom deterioration. This was done by sealing the end cap, 


which presses the target disc directly against the O-ring 


vacuum seal, with a plate into which short lengths of tubing 


had been brazed. (See Figure 3.3) Water forced through 


the cap removed the heat from the target disc. 


na 


peo ~=ClLhhe Neutron Source 


As in previous experimental measurements of neutron 
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Slows cova time (1,2), the Tae Cone reaction was used 
as the neutron source mainly on the Dasis of item en | 
yield and relatively good homogeneity. For a given.angle 
meomeune vargeu and for laboratory = provom energy UD 50 2. 34m 
Mey, une reaction is monoenergetic. The threshold of the 
Li’(p,n)Be’ “reaction is shown in Figure 4.6 of Adamantiades’ 
ee Because of neutron intensity considerations, a 
proton energy of 2.4 MeV was used in this experiment. 

The target used in these experiments consisted of a 
layer 5 Saseeoc al lithium 0.8/cm* thick deposited through 
evaporation an a Circular Gantvalumee tec 0.055 inches thick 
mgt 5/16 inch in diameter. Although the thickness of the 
target feeee eee he ene the range of 2.5 MeV protons 
re lithiun, some protons will undergo slowing down colli- 
Sions before undergoing nuclear reactions. The energy 
spread of the actual beam of protons incident on the tar- 
pee sowever, is very small in ie Rockefeller Van de Graaff 
pereclerator. 

| To eee: there were three major sources of neutron 
- energy inhomogeneity in this experiment: 

‘1), Angular dependence of neutron energy. 

2) Excited product nucleus reaction. 

3) Energy spread of the reacting protons. 

As mentioned in beet came this energy inhomogeneity 
could be tolerated since the slowing down time is relatively 


e 
"ai 


Mescensitive to the source energy. 





B. BinesVoderatine: Ascembia: 
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The moderating assembly (see Figure 3.4) consisted of 
reactor-grade graphite stringers, approximately 4.0 inches 
by 4.¢c5 inches and 4.0 feet long, stacked so that alternate 
>» layers were at right angles to each other. The overall 
dimensions of the assembly was 3.6 feeb by 4. CMe cuap., 

4.0 feet. The entire assembly was covered with aluminum 


foil to shield the monitoring and detection systems from 


me RAE. pickup coming from the pulsing system, although this 


pickup was never entirely eliminated. 
aia In a stringer located at the center of one face of 
_ the pile, a hole was bored so that the target, along with the 
meecooling water tubes and the cables of sume monitoring sys- 
tem, could be Onaerted into the center of the assembly. 

At a right angle to, and on the same horizontal plane 
as the beam tube was an open channel, 4.0 inches wide, 
Seeo inches high, and 19.3 ener lowe” which led from the 
ae monmone face of ene wescenbl It .was through this 
channel ae the detector, foil, and graphite spacers were 


inserted or removed (see Secti1 on. 


me Lectronic Systems 
joe =6vBeam Pulsing System 





Peepnouehn the Rockefeller accelerator nas yamvcrnainaL 
PUlSing capability, its single sweeping frequency was much 
moo high for this experiment. ‘therefore, a post-accelera— 


tion system was used which had the advantage of a variable 
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pulse width and repetition rate (see Figure 3.5). 4 
detailed discussion of the system is given in Reference (1). 
A repetition rate of 2kc/sec was used, allowing 500 usec 
between pulses, which was enough time to allow ale: tran- 
sient neutron levels to die out. A pulse width of 200 nsec 
was used in this experiment. 

The pulsing system was triggered by the General Radio 
meme cenerator (type no. 1237-A), which was alsouusccamuce 
start the timing circuits of the analyzing system (see 
feevion 4.5.3). The output pulse of the G.R. pulser 
triggered an 2-H Research Laboratories Model No. 130 pulser. 

The variable width output pulse of the E-H pulser was 
the input for the driver circuit which is essentially a 
pulse amplifier. An extensive discussion of We antes 
circuit is given by A.E. Waltar (3), who designed it. 

’ Vaughan modified the circuit somewhat to give a narrower 


pulse. in this experiment the driver circuit delivered a 


eae keV pulse to the deflection plates. 


feoee Che Miltoctienecveten 
The object of these experiments was to detect neutron 
capture gamma rays as a function of time, so an efficient 
gamma-ray detector was required. A Harshaw Integral Line 
Seoul liator Assembly, consisting of a 3 inch kK 43 inch 
Weal(Tl) crystal attached to a 6363 DuMont Shores 


tube, was used. The crystal and tube are hermetrically 


sealed within an aluminum housing, making the detector 
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a simple plug-in unit. «a thermal neutron shield, consisting 
of 4mm. thick boron carbide enclosed by aluminum, was pro- 
vided to keep activation of the sodium at a minimun. 
The solid state preamplifier, which plugged into the 
,end of the detector assembly, was designed ane spice 
cylindrical can having the same diameter as the detector. 
Vaughan gives a schematic diagram of the preamplifier in 
Figure 4.6 of his thesis. 
The high voltage for the phototube was provided by 
a 2 KV power supply, wnile the preamplifier was powered by 
asmall solid state power supply. The entire detector 
Was wrapped in-aluminum foil, and all leads to and from 
tne detector were enclosed within copper braid to eliminate 


J 


fee. pickup. 


mee = the Analyzine System a, J 

The main features of the time analyzing system were: 

1) A start pulse providing a time reference, derived 
from the beam pulsing. network. 

2) <A stop pulse sweneeare the time at which an event 
Seeemcerest took place, eee from the detector. 

4) A timing unit, used in conjunction with a multi- 
channel analyzer, whicn produced a count in ee enamel 
erresponding to the time interval between the arrival of 
the start and stop pulses. Conventional, commercially | 


available equirment was used in the entire analyzing system 


(see Figure 3.36). 
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pulse was obtained from the Gencral Radio 
Brigeer (see Section 3.5.1). This pulse Was fedunmcuce 
a General Radio 1392-4 Time Delay Generator which provided 
peotabdle, variable time delay. Unfortunavely  eciees ioe 
exise of the G.R. delay was too narrow Go Griecer Giese 
ing unit, and an integrating cae employed to 
stretch the vulse. 

The detector output signal was fed to two Ortec 410 
pm robbed Discriminator, with the strobe input coming from | 
an Ortec 407 Gero Cross-over Pickoff Unit. The cross- ! 
over pickoff was fed by the other Ortec 410, whose time con- 
stants had been optimized for a stable cross-over point. | 
A positive logic signal was provided by the discriminator | 


(} 
when the incoming detector signal met minimum energy require- 


| 
meee 6 ce strobe assured a constant time relation betweeu 


M@eecetbector input and discriminator output. The positive 


—- 


logic pulse was attenuated by a resistor network before it 
ae 
was fed into the timing unit. 


The timing unit consisted of a TMC Model 211 Time of 


aoe + 
os ee Se 


Flight Logic Unit which plugged into the TMC Model CN-110 
256 Channel Pulse Analyzer. The 2il unit converted the Anne 
ference between the trigger input and a detector-signal 
tO an appropriate channel number, depending on the scale 
/tecor which could be varied from 0.25 to 64 isec/channel. 
This information was then stored in the multi-channel analy- 
Berl. | | 


The TMC cll Unit could nandle more than one stop pulse 





(detector signal) per cycle. However, lope o ew 
mecded sbeo-analyze anil record e¢acn incoming sisnaley.s 
correction for this deadtime was included in the computer 
program used to process the data (see Section 4.3). More 


a 


detailed information on the 211 unit and the multi-channel 


analyzer is ies in Reference (4). 
me. t ine WMonitorine and Normalizing. s ys Lem 

To insure that the pulsing system was operating under 
the best possible conditions at all times a beam pulse 
Monitor was provided. 

The beam pulse monitor was necessary due to the pro- 
blem of focusing the proton beam from the accelerator. 
Because of various changes and dr aGS aimee ile accelerator, 
the setting for optimum focusing changed irom day stosaa 


and had to be constantly adjusted. In pulsed operation, 


tne only practical method of focusing was to use an electri- 


abe 


Cal signal from the target itself and to observe the signal 


Mith a fast oscilloscope. The system employed is sketched 


Moerigure 4.7. Notice that the target was insulated from 


the beam tube with a glass ring which enables the signal 


tO be picked off the target. The monitoring system was | 
Sesvecially susceptible to R.F. pick-up from the pulser. 


All cables of the system were shielded against pick-up; 


Mmencne problem was never completely eliminated. A typical 


bareet pulse can be seen in Figure a 1d of Vaughan’s thesis. 


The electronic pulse being 
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= 2 o meee : a gt pes Perey Mm BS = S . ; 
Could a@kso De monitored. This system is alsowsie7 


— 


BaeSUTre 9x '/s 
A meuGtron monitor was used as a means Cf merical. le 


° 


foe comparing different runs and subtracting Uae roma 
meaetailed discussion of normalization is Ziven tn Cnapees 
5 of Reference (1). As in previous experiments, this author 
concluded tnat toval neutrons produced during an experiment 
mould be the most reliable quantity to use of normallzationms 
A boron triflouride counter was placed between the graphite. 
stringers of the assembly about 13 inches from the target 
and on the side opposite the absorber foil and detector. 
mous, With the output ‘of the a tube going to a scaler, 
the time integrated neutron.flux as seen by the detector 

i} 
could be recordc?. 


? : 
' t 
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CHAPTER 4 


Experimental Results and Discussion 


4A.j) The Use of Resonance Absorbers as Neutron Density 
indi cavors 

Resonance absorbers were used in this experiment to 
sample meu ron ee the resonance energy. The 
most suitable materials for this purpose are those which 
exhibit a strong isolated absorption resonance. The 
nigher and narrower the resonance, the better it eee 
mates 2 delta function. Table II eG Reference (1) Fives 
a list of absorbers which exhibit strong isolated reson- 
meances. The energy and peak cross Section of Stieeme oe. 
Meee is peso Given. 

mpoovner important consideration, in ehoosimg sa .cuesr—as 
able ee ber was its prompt neutron capture gamma-ray 
spectrun. This did not prove to be a problem since the 
de-excitation gamma rays of most materials are spread 
over a wide range of energy. Reference (1) gives capture 
gamma-ray spectra for the three absorbers used in these 
experiments: cadmium, indium, and gold. The capture gamma 
Spectra of these materials extends from the KeV region = , 
to ee exceeding 7 MeV. | 


Hndium and gold were chosen for their strong asolated 


resonances at 1.46 eV and 5.0 eV respectively. Cadmiun, 








vnich has a broad resonance that reaches agmaiauiiae 
Maabout O.c eV, was also used, although it was expected that 
the time resolution for these measurements would not be 


me 2O0GeeSs vRauy tor:-indium andycoid: 


im order to insure @ constant @eometricure aviononume. 
Meeveocn the foil and detector, a foil holder was used. 
The foil holder consisted of a one inch thick slab of era- 
phite, which was sized to fit into’ the channel of the assem- 
Bly and was bonded to the end of the detector with epoxy 
cement. The foil was taped to this slab, and the whole 
Mebecvor-foil array was inserted into the graphite assen-— 
bly through the channel (see Figure 4.1). \ | 

me To measure the spatial dependence of the slowing down =: 
Memon it was necessary to be able to move the detector— 
foil array along the cnannel. Graphite slabs of varying ) 
thickness, which were inserted in the channel between the | 
source and foil, assured that there would always be solid 3 
M@eoeraling material in front of the foil. | 

Beiore performing the actual €xperiments, short: runs. +. | 
were made while varying the discriminator of the Ortec : 
413 Strobbed Discriminator. For each foil, the discrimin- ~ : 
ator setting which resulted in the best foil counts to : 
background .counts ratio was determined. | 

| A typical set of exverimental data is shown in Figure 
\ ; 


oe a)a, These results are for a 7 cm x 7 em we. e> cr gold ’ 


2 } . 
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foil positioned 23 cm away from the target. The bdeam 

- . | 
pulse width was 200 nsec, and the time interval between 
bursts was 500 ysec. The channel width was 250 nsec. 
Segime distribution wi h the s0ld absorber” ing placem i a. 

- taken, and a distribution witnout the adsorber was taken 
immedia&ely afterwards. The Bi counts for the goldeiows 
was 251,111 and for the background was 259,992. Notice 
the very prominent peak around channel 20. This peak is 
called the prompt peak and occurs when the burst of pro- 
mons hits the lithium target. This peak was used as an 

saccurate zero time reference in the experiment. The right 
Mart Of the sraph (notice the change of scale) shows the 
region in which the most probable slowing down time occurred. 


' 
No peak is discernible from the raw data. 


4, 


keround Subtraction, Dead Time 





, Pemmentioned in Section 9.5.4, the Gime integracved 
m@emmoron tlux, as derived from a BP, tube, WaS USeCd fiom 
normalization. The number of BPs counts registered by 
the scaler was recorded for each experiment. A factor 
obtained by taking the ratio of background-run BPs counts 
co Foil-run 3F, counts was applied to the foil data. 
| in order to obtain the net effect of neutron cap- 
Mees in tne foil, the counts originating from all other 


cod 


PeyenusS MUSt be sud ec when the time divstrereucion 24 


taken. These sources of parasitic gamma rays include: 
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1) Inelastic scattering by the target, the absorber 
Polen Foe de vector. 

o ) - Absorption in materials other than thesabecrpe: 
foil. | 


¢ 


i %) Reactions coincident with tne neutron onec 


YR 


wy y 


4) Overlapping events from preceding bursts. 

B) steady state sackeround. 

Ali these items were taxen into account by taking 
a time distribution witnout the absorber foil in place, 
and then subtracting 1t from the distribution taken with 
Mae foil. 

As mentioned in Bocce 5.5.5, there was a 16 usec 
dead time associated with each count registered in the 

(3 

Multi-channel analyzer. In Appendix C of Reference (2), 


a theor: 


\ 


tical discussion of the dead time correction for 
the TMC Model 211 Unit is presented. The correction fac- 
bOrs derived in that reference were used in this thesis. 

Computer programs Wile nostenneees subtracted the 
background, and applied dead time corrections to the raw 
data are listed in Appendix A. 

Pigure 4.e(b) shows the time distribution resulting 
from processing the raw data in Figure 4.2(a). It is 
obvious from this graph that there is a maximum in the 
Mime cistribution occurring roughiy at channel number 60, 
[e's also obvious, however, that the fluctuationseia the 
curve prevent an accurate determination of the maximum 


Moincg., lo reduce this uncertainty, a computer program 
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developed at #H.1.7. (5) which applied a Smcotmimene oe 


\ 


mique to the corrected data was used. ~Hifecvives  aamee 


morogran fits a fouler series curve Uovtne  daldeporan 

; | 
mpmoothine 1S accomplished by CchoosSin@ not Go Tit suiemeea. 
gPigh frequency waves, which correspond here to sta- 


Histis 41 fluctuations. The frequency .cutoff, above which 


all big er frequencies will be thrown aay ep toe nee ame 






contr iM of the programmer. This ean, cesieeed sdearoc mem 
smoottne can be effected Go a seu of davgavel mcs ee 
bees experiment, all frequencies less than four channels 
Beege were drooped. Figure 4.2(c) is the smoothed version 
Si figure 4.2(>). Using this smoothed distribution; it 
was estimated that the maximum point of the time dis- 


tribution occurred in channel Camo beet ec the prompt 
Mea. is in channel 20, the slowing down time is: (43.5+1) @ 
x (.25 psec/channel) = 10.87+.25 sec. The smoothing 


“program is listed in Appendix A. 


4.4 Snace and BEnerev Devnendent Slowing Down Time Mea- 








Surements 

If various absorbers with resonances located at 
different energies are used and if the distance of these 
meoorbers irom the source is varied, the time dependent 
neutron distribution can be measured as a function of 
space -oomm@euceon enerk 7. 


Fieures 4.3 to 4.6 represent time distributions to 


the cadmium resonance, Figures 4.7 to 4.11 represent time 





i 


(qq) @°t aanstyq Jo uotsro, paygooug (9) 2° aansty 


| BABWAN TaNNGHS , . 
Ole D2 N12 coi oss Az 1 - 096 oe 7 


Cus 


2 
Le 


LO) 


—ree. 





‘ ~“ ~~ 
cd 


= 
a 


WOLI 





~6- 


distribution to the indium resonance; and Pisurecm ae 


‘to 4.16 represent time distributions to the gold reson- 


- 


ance. In each Figfure, (a) represents the distribution 


merore smoothing and (b) after smoothing. Notice than 


'the scale for the cadmium curves is 500 nsec/channel, 


while it is 250 nsec/channel for the indium and gold 

The results of this e€Xperiment, vosecvner Wives e mc im 
probable Crees, Powe oie? presented in, Figures 
4.17 to 4.19. The experimental uncertainty associated 
feeeecach point reflects the broadness in the time dis-— 
tribution: around the maximum value. As expected, the 
peaks of the cadmium distribution were much broader than 
Ghose of gold or indium. The solid line in the Pomel and 
moe figures has been constructed from theoretical time 7@ 


Mester ibutions (see section 2.4), based on diffusion theory 


(4). 


The experimental results for.the slowing down time 
jem.) eV appear to be in rood agreement with the theo- 
retical calculations. The agreement with theory for 
1.46 eV is not as good. The experimental results indi- 
Mee onat there is little, if any, spatial dependence of 
the Slowing down time for distances closer than about 
40 cm from the source. For distances greater tnan 30 cm, 
the experimental results show a snarp rise in tes witn dis- 
memee, indicating a strong spatial dependence.  Theore- 


oe calculations do not PredliCumas Silatp eae ce. 
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theory and experiment was 


ce 
Vv 


probably due to the breakdown of the infinite medium 
Boundary conditions, since the foil was only 9 cm from 
the outer edge of the assembly at the last spatial point. 
The absence of moderating material surrounding Che wa ous 
meant that the neutrons took a long2r time to slow down 


a 


to a given energy. More measurements need to be taken 


mumecreaver distances to resolve this question. 


A few measurenents were made witn the detector separ- 
aoed from the toil by abdout 10 cm of graphite in order 
momaetermine if the perturbation caused by the detector 
Beoesany effect on the slowing down time. Theselresults 
are marked by triangles in the Figures. It can - Seen 
that the gold and cadmium results agree, within the limits 
meee DeCrimenvcal uncertaintity, with previous measurements. 
The indium point, however, falls far below the original 
measurements. The author can give no explanation for the 
ee adiction of this point with eiuher the crigiam results 
for indium or the close agreement between the ee and 
second measurements for the gold and cadmiun. veo 


Meesirements, with the detector and foil separated, need 


mombe performed to clear up this contradiction. 


fee COnclusions 


There was some agreement between the spatially depen- 


moo Slowing down time measurements made an Gnese experi 
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Mencs and theoretical calculations using diftiusion theory, 
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No measurable spatial dependence of the slowing down time 


was observed for distances less than 30 cm. Spatial depen- 


-77- 


dence was observed for cistances greater than 50 Cia em 


this dependence was much more promounced in the experiments 


than predicted by slowing down theory. More space depen- 


a 


ent measurements at distances greater than 44 cm should 
be made in a larger sraphite assembly to determine the 
Mkiemopatial dependence. Measuremenvs to determine ae 


Mee cruuroation of tne detector nad any eliecco on tie 


feowins down time were inconclusive and should be repeateds 
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Appendix A. Computer Program Listings 


MORMALIZATION, DEAD TIME CORRECTION, AND BACKGROUND 
SUB ERACTLON CODE inbST ING 


DIMENSION T8M(255) ¢CHAN( 25S) »BACK( 255) »FOIL(255),ANORM( 255), 
LCBACK(235) ,CFOILI255),X(255) 9¥(255)s3TBA(236),ACHAN( 206) ,B{ 800) 
COMMON TBM,CHAN,WIOTH ! 
FOULVALENCE(193M(50),TBACLI 29 (CHAN(50) ,ACHAN(1)) 
READ 200,WIDIH 
200 FORMAT{F5.2) 
READ 99,8COUNT,FCOUNT | 
99 FORMAT (2F10.0) — j 
NUM=0 | i 
103 NUM=NUM#lL 
- READ 100;(BACK(1),1=1,255;3 
100 FORMAT (3X;10F7.0) 
READ LOL, (FOIL(1)41=1,255) 
LOL FORMAT {3X;10F7.0) 
READ 98,88UR,ABUR 
98 FORMAT (2F160.0) | 
PRINT 7O0l,BCOUNT,FCOUNT U 
TOL FORMAT (3X%,2F1G.0) 
PRINT 100,{BACK{(1),1=1,255) 
PRINT 101,{FOIL{!),£=1,259) 
. PRINT 702,8BUR8,ABUR 
Ho? FORMAT (3X72F10.0) 
L=-] 
ea20,L=L+1 
he £19) 419,419,417 
419 00 418 1=1-+255 
418 X(1)=BACK(1) 
BUR=BBUR 7 
GON TO 415 7 
417 DO 416 1=1:;255 
416 X{I)=FOIL(1) 
RBUR=ABUR 
415 DUM=16./WIDTH 2 
M=DUM ; 
Y(L)=BURFLOGFIBUR/(BUR-X{1))) , : 
SUM=0. 
DO 610 T=2,N 
SUM=SUMt#X{(IT-L) | | 


mm eee eee 





S = Pgs 


OOM 


610 


615 
414 
413 


412 
411 


906 


306 


30 


T00 


Y( 1) =BURsLOGF((BUR-SUM)/(BUR-SUM-X{(I))) 
SUM=SUN#X{N) 
NM=\+ 1 
OO 685 [=NMy255 
Recs nes cauR-Sus 1s thee 
K={i-N 
SUM=SUM4X01T)-X(K) 
IF (b)4147414;412 
DO 413 T=1,255 
More K LL )=¥(1} 
GO TQ 420 
DO 411 L=1,255 
@poti{1i)=y{tt} 
Z=BCOUNT /FCOUNT 
00 906 [=1,255 
Maoem( il }=ZaCFOIL{( i} 
TBM( TP =ANORM(CT)-CBACK(T} 
BPANtTI}=FLOATF(I) 
CONTINUE 


PRINT 306,(TBM(I),10T=1,255) 


FORMAT [1H ,LOFL0.1) 
CALCOMP ROUTINES 


Batt PLOTS1(8,800) 


BALL Wes a eee ID.NOw = M3994-5788,0.,20) 


PALL MS21(4.078.5,0¢2sNUMe 00 y-1} 
mee ALL ate. 5274e2,15H CHANNEL NUMBER, 


pO lPILeE 10 


READ 99,;BCOUNT,FCOUNT 
IF (BCOUNT) 70077004103 
Pee EX{LT 

END 


foes i 9H COUNTS PER a 19,ACHAN, Es 206,0-30) 





mon 


200 


102 


201 


205 
203 
99 
103 
100 
101 
9g 


TOL 


as). 2 
420 
419 
418 


41.7 
416 


4 


DATA SMOOTHING CODE LISTING INCLUDING NORMALIZATION, DEAD 
TIME CORRECTION, AND BACKGROUND SUBYRACTION 


DIMENSION TBM{37)sXSIN{75) sSEIGEN(373 75) ,KXCOS1 (15) 4 %GO5207 oe 
LXCOS4(37},78M(255)eTBK (255) »~CHAN(255).THETA(75}),ALPHA(75) » 
ZBETA(75),GAMNA(75),W7{75)-B(800),BACK(255),FOIL(255) ,ANORM(255),y 
BCBACK(255) ,CFOIL(255) »X(255) 5 Y(2595 II BAL 21 Ong AOrAN oclon 

COMMON TBN,XSIN,SEIGEN?;XCOS1,XCOS2,XCOS4,1TBM,TBKsCHAN;, NIDTH 

EQUIVALE WEE TBR Oe BT ae 

READ 200;,IW,OMEGAH,CRI,FACTR WIDTH 7 

FORMAT (14,610.49 F60.93F6.3256-2) ee 

Meee O2,(W1( 1} ,l=lstiwW) 

FORMAT (12F6.3) 

PRINT 2015 0W,OMEGAH.CRT;,FACTR,WIDTH 

FORMAT {(LHl-SH N-OMEGA=,14,1X%99H D-OMEGA=,E10.4,1X, 

PS CRT=,°6.0; i SH FACTOR=,F6.3,7H WIDTH=sF6.2) 

PRINT 205 ~ 

FORMAT (1H ,»20X,16H WEIGHT FUNCTION) a 

PRINT 203,(WT(1},1=1,1W)} es 

FORMAT (LH »5615.7) : , 

READ 99,BCOUNT,FCOUNT 

FORMAT {2Fi0.0) 

NUM =O 

NUM=NUM+41. ~~ ne ot 

PEMD 100:(B8ACK(1),1=1,255)} . 

FORMAT (3X,1iGF7.0} . | | ° 

moro 101, {FUIL(I},1l=].255)} 

FORMAT (3%s10F7.0}) 

READ ¥8,B8UR,ABUR 

FORMAT (2F10.0} 

PRINT 7FO1;3COUNT,FCOUNT 

FORMAT (3X,2F10.0) 

Pein t 100, (MATKII) «Lei anZ255} 

PRINT LO1,(FOIL{(L),i=1,255)} 

PRINT 702,B8UR,ABUR 

FORMAT (3X,2F10.90) 

b= — | , ; 2 

t=Lel : 

Pe) 111419,419,417 

Pemels {=1,255 

X(I)=BACK{(1) 

BUR =BBUR 

GO TO 415 

oeera. [=],255 

MpierOiL{t{{) 

BUR=ABUR 





415 OUMN=L6./WIDTH 
N=DUN 
Y (1) =BUR*LOGF{BUR/{BUR-X(1})) 
Bens 
00 610 I1=2,N 
SUM=SUM4tXC1-L) 
6L0 Y(i}=BUS*LOGE( (BUR-SUM) /{BUR-SUM-X(I)))} , 
SUM=SUMEXIN) | : 
NMoN 4] 
00 615 I=NM,255 
¥(1)=BUR#LOGF((BUR- SUKI / (BURG SUM~X{1I))) | 
K=I-N | 
Bi5 SUM=SUM+X(1)-XtK) | ! 
Pe 11.9414541454512 
4144 00 413 1=1,255 
fs CBACK(1I)=Y(1) 
60 10 420 
412 00 411 I=1,255 
MipecFO{L(1)=Ytl) 
DO 905 I=1;37 i" 
YKAKUG=OMEGAHSFLOATF(IW-1) #FLOATF{(1) 7 ' 
XCOS4(1)=COSF(YKAKU4) : 
905 CONTINUE 
Bo 9Ol 1=1,37 | 
00 902 J=1l,IW | 
XKAKUSOMEGAH*FLOATE( IL) #FLOATFE( J-1) 
eEVGEN (1, JJ=SINF(XKAKU) 
902 CONTINUE 
901 CONTINUE ! 
oo 90% Ja1,1W° : 
THETA( J) =OMEGAH#*®FELOATF(J-1) 
ALPHA( JP=THETALJS) ¥#(-3) #(THETA( JS) *##2¢THETA( J) #SEIGEN(L, J). 
LCOSFI(THETA(JS))~-2. #SEIGEN(1,J)##2) 
BETA(J)=20#THETA(I]##(-3)"(THETA(J) #lLe#COSE{THETA(S) }*#2)-2e% 
LSESGEN{1L+J) eCOSFE(THETA(J))) 
GAMMAL 3) =4.2THETAL JS} #*(-3) #(SEIGEN(12J)—- THETAC J) ®COSFITHETA( JS) 1) 
903 CONTINUE | . 
DQ 906 1=1,255 i 
Z=B8COUNT/FCOUNT | — 7 
ANORM{LY=Z2CFOILIL) ° 
TBM{ 1) =ANORM(L)-CBACK{T} | 
906 CONTINUE | . 
00 15 I=1;27 | 


—_— = = — om . _ 


| 
y 





511 


517 


Oe 904 J=)41 

YKAKUL=OMEGAH#FLOATF(J-1) 

PAR I2eZ TSG YRAKUL 

XCOSL (JS) =COSFLYKAKUL ) 

mCOSZ(J)=COSFE(YKAKUZ2} 

Pot i NUE 

09 511 I=1,1f# 

YA=ALPHA(T) 

Pee hepAt{ i) 

Pe=GAANAL 1) 

S$1G64=0. 

SIG3=9,5F(TBN{(LI#SEIGEN{1, T}+TBN(27)#SEIGEN( 27, T}) 

DO Ll NN=3,25;,2 

SLG3=SIG3¢+TBN(NN) SSEIGEN(NN, 1) 

Ba 2 NW=2,26,2 

SIG4=SIG4+TBN(NN) *SEITGEN(NN,T) © 

RAH=VYAR( TENG LISXCOSLETI-TEN(27) *XCOS201))+Y8eSIG3+YC2S1G4 
ASINCLI=XA#y7T (1) | 

CONTINUE \ 

Pie Si? [X=1,16 
MA=ALPHACIX+#1L) 
YB=BETACIX+]) 
YC=GAHMALT XK+1) 
S$1G4=0. 
SIG3=O0.52{({XSINCLISSEIGEN(CIA,LIFXSIN( IW) FSELCGEN(IX,IW)) 
NY={L4-2 


0a 3 NN=3;NY;72 


SIG3=SIG3+XSIN€NN) ©SEIGENCLEX,NN) 

NY=I[W-1 

DO 4 NN=2,NY-2 

SIG4=SIG4+XSIN(NM) #SEIGENLIXyNN) 
XA=OMEGAH#(YA®(XSIN(LI-XSIN( IW) *XCOS4( 1X) )+YBaeSIG3+YC#S1G4)} 
IN=IX 

TBK{ IN) =0.636622K%A 


 CHANCIN}=FLOATFCIN) 


CONTINUE 

DO 914 J=1.TH 

YKAKUL=OMEGAH#FLOATFLJ-1) | | 
¥KAKU2=37.2#YKAKUL . 
ACOSL( JS) =COSFLYKAKUL) | 
XCOS2{J)=COSFLYKAKU2) 





525 


ac 


527 
530 


935 


eee 


CONTINUE 

00 530 J=2,;,15 

DO 525 1=hs27 

K=l6e(J-1)-LOF] 

TEN(T)=TS8M(K) 

CQ $21 =1L,1W 

mA ALPHAI I } 

YB=ARETALL) 

YC=GANMA(T1) 

mG 4=0. 
S(63=0.52(TBN(L)I®SELGEN(L, 1)4+TB8N(37) sSEIGEN(3731)) 
D0 5 NN=3-35,2 
SIG3=SIG3+TBN(NN} ©SEIGEN(NNG TL) 

DO 6 NN=2,36;72 

SiG42SI1G4+T3NONN) ®SEIGEN(NNeI} 

XA= YAR {TBN(L) #XCOSL(I)-TBN(37) #XCOS2(1))#YB#SIG3+YC#SIG4 
MSIN(LI=XARaT(L) 

CONTINUE © 

M527 1X=11,26 

YA=ALPHA(IX+1} | \ 
YB=BETA( I1%+1) 

YC=GAMMA{IX41} 

51G4=0. 

1G3=0. S2{XSIN(L) #SEIGEN(IX, LI +XSINCEW)#SEIGEN(IX, TW)) 
NY=twW-2 

DO T NN=3,NY,2 | 

SIG3=SIG3+XSIN(NN) #©SEIGEN(IXENN) | 

NY= tW~-] 

DO 8 NN=2,NY22 

SIG4=SIG4+4XSIN(NN} #SELGENCIXyNN) 

XA=ONEGAH* (Y hs (XSTNCL)“XSINCIW) FXCOS4( 1X1) #YB*SIG3+YC#S164) 
[N=162#(J-L)}+1X-10 

TBXLIN) =0.6365624XA 


CHAN(IN) SFLOATFCIN) 


CONTINUE 

CONTINUE 

OO 535 1=1,25 

K=230+] ‘ 
TBN(L)=TBM(K) 3 

Be 924 J=1,I1n 

YKAKULSOMEGAH2ELOATE(J-L) 
YYAKUZ=25. 8" YKAYUL 

XCOSL{ J) =COSFLYKAKUL) 

XCOS2( 4) =COSELYKAKU2} 





MmOoOnm 


92% 


10 


3 1 


11 


2 


237 
307 


306 


2Q 


CONTINUE 

DO 53) [=lL,10W 
AP EA tT} 
YS=BETACT1} 
Re=GAMaA i) 


rie = 0) 
5 wey SS a 


SIG3=O0.5#(TBN(LI®SEIGEN(1, £)+TBN(25) eSELGEN(25,1)) 


DO 39 NN=3;2372 
STG3=SIG3tTBN(NN) *«SEITGEN(NN,T) 
DG 10 NN=2424%%2 
SIG4=S1G64+TBN(NN) #SEITGEN(NN. 1) 


XA=YAPLTBNGILISXCOSL CIT) -TBN(25) #XCOS20CT} VtYBaSiG3+YCeS1G4 


ASINCT) =XA#enwT (2) 

CONTINUE 

. eof ixX=1lisy25 
ALPHA IX+]) 

OE ast 3 

YC=GAMMACIX+#1) 

$164=0. 


Ca 


91G3=0, De (XSIN(L) ®SETGEM GLX) 1) #XSIN(TW) SETGEN( 1X4 1H1) 


NY=iW-2 

DO Lh NN=35NY22 

SIG3= SEGS+XSENCTHNI #SELGENECTXSNNG 
Ny=[¥— ]. 

DO LA NN=2eNY22 

SEG =SiG4+XSTN(NN) #SETGENCEXyNN)} 


| 
XASOMEGAH® (YAR (XSINCLI-XSINC IW) ®&XCOS4( IX) ) +YB2SIG3+YCeSIG4) 


IN=Z2290+14 

TBKCIN)=0.63652#KA 
CHAN(IN)J=FLOATFCIN) 

CONTINUE 

PRINT 307 

BIRMNAT {1H1,18H SMODTHED FUNGTION: 
PRENT 306, (78X{1)},T=1,255) 

FORMAT {1H ;1LOFLO.1L) 


CALCOMP ROUTINES 
CAE PLOTSL(B;3001 


CALL NUMBR1L64.0+8.57002sNUM,06.9-1) 
GALL S7TU8BL5 (3.06930 vere ID NO. 


M3994-5788,0.720) 


CALL PICTUR(6.57%22;15H CHANNEL NUMBER, 


eOeriLle 10 


READ 99, SCOUNT, ECOUNT 
[TF (8COUNT) 700,700,103 


115.28H SMOOTHED COUNTS PER CHANNEL »28,ACHAN,TBA, oe 





OO CALL EXIT 


END 
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